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THRESHOLDS FOR REDUCTION-RELATED ENTANGLEMENT CRITERIA IN 

QUANTUM INFORMATION THEORY 

MARIA ANASTASIA JIVULESCU, NICOLAE LUPA, AND ION NECHITA 


Abstract. We consider random bipartite quantum states obtained by tracing out one subsystem from a 
random, uniformly distributed, tripartite pure quantum state. We compute thresholds for the dimension of 
the system being traced out, so that the resulting bipartite quantum state satisfies the reduction criterion 
in different asymptotic regimes. We consider as well the basis-independent version of the reduction criterion 
(the absolute reduction criterion), computing thresholds for the corresponding eigenvalue sets. We do the 
same for other sets relevant in the study of absolute separability, using techniques from random matrix 
theory. Finally, we gather and compare the known values for the thresholds corresponding to different 
entanglement criteria, and conclude with a list of open questions. 
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1. Introduction 

The notion of quantum entanglement has been proved to be at the core of many quantum phenomena, 
such as teleportation, dense coding or cryptography. Moreover, it is a key ingredient to the computational 
power of quantum devices. Entanglement expresses inseparability, that is unusual correlations between the 
subsystems of a quantum system which cannot be explained by classical, non-quantum, models. Hence, 
a central question in the theory of quantum information and computation is detecting the presence and 
measuring the amount of entanglement present in a given quantum system. 

One of the most efhcient tools in detecting entanglement is the positive partial transpose (PPT) criterion 
[26]. It states that if a quantum state is separable, then the partial transpose with respect to one of the 
subsystems is positive-semidefinite. It represents a necessary condition for separability and more often it is 
applied as a tool to detect entanglement: if the partial transpose of a given state is not positive-semidefinite, 
then the state is entangled. As the partial transposition criterion is obtained by applying the transposition 
operator over the second subsystem, it raised the question of finding other positive maps P with the property 
(id(g)P)(p) is not positive-semidefinite, for some entangled bipartite states p. One needs to look for such maps 
P in the class of positive, but not completely positive applications. In [14], the authors show the following 
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converse: a quantum state p G A^„(C) 0 A^fc(C) is separable if and only if (id(g)P)(p) is positive-semidefinite 
for all positive maps P : Mk{C) —>■ Mm{C) and all positive integers m £ N (actually, m = n suffices, see 
[14, Theorem 2]), where id is the identity map. Thus, each fixed positive map yields a necessary condition 
for separability; in other words, each fixed positive map yields an entanglement criterion. Note that these 
conditions (resp. criteria) become trivial for completely positive maps. 

A possible choice of the positive map P is the reduction map R : Mk{C) —t A4fe(C), R{X) := /fe-Tr(A)—A, 
and the corresponding separability test is called reduction (RED) criterion [9, 16]. The reduction criterion 
is weaker than the PPT criterion: if a state violates the reduction criterion, then it also violates the PPT 
criterion [16]. Conversely, there exist states (some entangled Werner states [32]) which satisfy the reduction 
criterion but violate the PPT criterion. The two criteria are equivalent if the subsystem on which the 
reduction map is applied is a qubit [9]. 

Both PPT and reduction criteria are efficient theoretical tools in detecting quantum entanglement, al¬ 
though the former presents the obvious advantage of a more elegant form, which in addition requires less 
computations. In practice however, since impure entanglement is produced, the concept of distillation was 
introduced as the process to produce a pure maximally entangled state by local quantum operations and 
classical communication, from many copies of an arbitrary entangled state (see [7], [17, pp. 870]). Horodecki 
proved that a PPT state is necessarily undistillable [15]. This result sheds light on the fact that in high 
dimensions there are entangled states which cannot be distilled. These states, namely PPT entangled states, 
are called bound entangled states, contrary to entangled states which can be distilled. It is possible to show 
that reduction criterion and entanglement distillation are connected: any state which violates the reduction 
criterion is distillable; conversely, if a state can be distilled by a certain protocol, then the state violates 
the reduction criterion [16]. This result justifies the use of the reduction criterion, even if, from a purely 
entanglement-detection perspective, it is weaker than the PPT criterion. 

The separability problem was also approached by studying the class of absolutely separable states (ASEP), 
i.e. states that remain separable under any global unitary transformation [22], that means to find conditions 
on the spectrum that characterize absolutely separable states (constraints on the eigenvalues of a state p 
guaranteeing that p is separable with respect to any decomposition of the corresponding product tensor space 
[21]). This problem was first fully solved in the qubit-qubit case in [29], and then in the qubit-qudit case in 
[20]. Furthermore, it is known that there is an Euclidean ball of known radius centered at the maximally- 
mixed state ® Ik) such that every state within this ball is separable [12] (see also [31]), meaning 

that any state within this ball is actually absolutely separable. However, there exist absolutely separable 
states outside of this ball [30, Appendix B]. In analogy to absolutely separable states, states which remain 
PPT/RED under any global unitary transformation are called absolutely PPT states (APPT) /absolutely 
RED states (ARED) [31]. Necessary and sufficient conditions on the spectrum of APPT-states are given 
in [13], in the form of a finite set (albeit exponentially large in the dimension) of linear matrix inequalities. 
For the case of ARED-states, necessary and sufficient conditions are given in the form of a infinite family of 
linear inequalities, which the spectrum has to verify [19]. 

In this paper we approach the problem of separability and absolute-separability from a different perspec¬ 
tive. We aim to derive thresholds for the reduction and absolute reduction criteria and to give a complete 
picture of threshold points for the class of entanglement criteria. The threshold point is defined in the 
following sense: given a random mixed state pab € Mni'C) 0 Alfe(C), obtained by partial tracing over C'* a 
uniformly distributed, pure quantum state x G C"(8)C^(8)C®, where the s-dimensional space is treated like an 
inaccessible environment, we ask for the probability that the state satisfies an entanglement criterion. When 
one (or both) of the system dimensions n and k are large, a threshold phenomenon occurs: if s ^ c • f{n, k), 
for some constant c > 0, or s is fixed, then there is a threshold value cq of the scaling parameter, such that 
the following holds: 

(1) for all c < Co, as dimension nk grows, the probability that pab satisfies the entangled criterion 
vanishes; 

(2) for all c > Co, as dimension nk grows, the probability that pab satisfies the entangled criterion 
converges to one. 

The threshold phenomenon was introduced by Aubrun to study the PPT criterion [2] . Our main contribution 
presented in this paper is to complete the computation of the thresholds for the reduction criterion given in 
[18] and to derive the threshold for the absolute reduction criterion, in different asymptotic regimes. 
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The paper is organized as follows: Sections 2 and 3 aim to introduce the main concepts and notations used 
in the paper. In Sections 4 and 5 we compute explicitly the value of thresholds for reduction and absolutely 
reduction criteria, in different asymptotic regimes. In Sections 6, 7, and 8 we derive thresholds for some sets 
that express certain conditions on probability vectors and which approximate the set of separable states. In 
the last section, we gather all the results about thresholds for different sets, and we present open questions 
related to this subject. 
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IN’s research has been supported by a von Humboldt fellowship and by the ANR projects OSQPI 2011 BSOl 
008 01, RMTQIT ANR-12-IS01-0001-01, and STOQ ANR-14-CF25-0003. 

2. Entanglement criteria 

In this paper the set of density operators (positive-semidefinite matrices of unit trace) acting on is 
denoted by Dd and for bipartite quantum systems on tensor product Hilbert space C"(8)C^ = we identify 
Dd with Dn^k, subscripts indicating the bipartition (n will denote the Hilbert space dimension of the first 
tensor factor and k that of the second one, and both n,k > 2). 

A density operator p S A4„(C) ® Mk{C) (here Mni'C) denotes the space of all n x n complex matrices) 
is called separable [32] if it can be written as 

P = YlP ^^ i^i ® fifi 

i 

for Pi > 0, J2iPi — 1) ^od for unit vectors G C", fi G (throughout the paper we will identify quantum 
states with their density matrices). The set of separable states [32, 17] in Z3„,/c is denoted by 

SEP„,fe := {p G Dn^k I P separable}. 

Efficient methods for explicit characterizations of SEP„ ^ are not known and for this reason upper and lower 
approximations are of interest [17]. 

On any matrix algebra A4d(C), we define the reduction map, 

R : Md{C) -G Md(C), R{X) := Id ■ Tr(A) - A, 

where Id denotes an identity matrix of the appropriate dimension (here, d) and Tr is the usual, unnormalized, 
matrix trace. From the definition, it follows that the map R is positive, i.e. R{X) > 0 whenever A > 0. 

For a bipartite matrix A = Xab G A4„(C) O A4fc(C) = Mnk{C), its reduction over the second subsystem 
(H) is denoted by 

:= (idO i?)(A^B) =Xa® Ik- Xab, 

where Xa ■= (id O Tr)(A) denotes the partial trace over (B) of the operator A = Xab- We write the 
transposition map on any matrix algebra Md{C) as 0, and we also write 0(A) = A^; we denote the partial 
transposition of a bipartite matrix A = Xab by 

A^^ := (id(g)0)(A). 

The composition of 0 with the completely positive map RQ : X i-G Id ■ Tr(A) — 0(A) is the reduction map 
R defined above; one says that the reduction map R is completely co-positive. 

Every positive map P on M.k{C) defines an entanglement criterion [14, 17]: if, for p G D^^k, the matrix 
(id 0 P){p) is not positive-semidefinite, then p is entangled. Specializing to the reduction map P = R, this 
becomes the reduction criterion [16, 9], which is also related to the distillability of the state in question 
[17]. Every bipartite state whose entanglement is detected by the reduction criterion is also detected by the 
partial transposition criterion [26, 14], which is the above criterion for the map P — Q-, this follows from the 
above mentioned representation of R as the composition of 0 with a completely positive map. 

The set of density operators p G Dn,k having positive reductions with respect to the second tensor factor 
for the fixed tensor decomposition Mnk{C) = Mn{C) 0 Mk{C) is denoted by 

RED„,fc := {p G D^^k I > 0}. 
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The entanglement criterion based on positive maps [14] implies the inclusion SEP„_fc C RED„ fc [16, 9]. 
Recall also that the set of states with positive partial transpose is 

PPT„,fc := {pe Dn^k\p^ > 0}. 

Note that, when k = 2, the reduction and the PPT criterion are equivalent [16, 9, 18], i.e. they detect 
entanglement for the same states, so that PPT „_2 = R.ED„ 2 ; in general, the reduction criterion is weaker: 
SEP„,fc C PPT„,fc C RED„,fe. Furthermore, it is well known that SEP„_fe = PPT„_fe whenever nk < 6 [14], 
Occasionally we will write RED instead of RED„ ^ etc., as the dimensions of the subsystems will be clear 
from the context most of the time. For a sketch of the different sets corresponding to the criteria described 
above and their inclusions, see Figure 1; the figure on the left contains the set RLN of states satisfying the 
realignment criterion [10] (also known as the computable cross-norm criterion [28]). 




Figure 1. On the left, subsets of ^ corresponding to entanglement criteria and their 
inclusions. On the right, subsets of corresponding to the “absolute” versions of the 
entanglement criteria, as well as the sets GER and SEPBALL, see Sections 7 and 8. 


Let us now introduce the “absolute” versions of the entanglement criteria above. We denote by Unk 
the set of unitary operators acting on C”^'. The set of states which remain RED under any global unitary 
transformation U G Unk is denoted by ARED absolutely RED”): 

ARED„,fe := {p G Dn,k\yU G Unk ■ {U pWf^'^ > 0} = f| [/RED„,feC/A 

ueu„k 

Similarly: 

APPT„,fe := f| UFPTn,kU*, 

U(iU„k 

ASEP„,fe := Pi USFFn,kU*. 

Obviously, ASEP„ fc C APPTn fc C ARED„ fc and AX„ fc C A„ fc, for X = SEP, PPT and RED. The question 
whether a quantum state p belongs in one of the three sets introduced above depends only on the spectrum 
of p; this is why sometimes we identify the sets AX with sets of spectra: 

nk 

AX„,fc C A„fc := {A GRf\Y,>^^ = !}• 

There are known results on the characterizations of the sets APPT„^fc and ARED„_fc as given by necessary 
and sufficient conditions in the form of families of linear inequalities which the spectrum has to verify [13, 19]. 
To date, there is no simple characterization of ASEP; in [20], Johnston shows that ASEP „_2 = APPTji _2 (and 
thus, also equal to ARED„^ 2 , see [19, Proposition 5.1]), while in [1] further evidence towards the conjecture 
that ASEP„_fc = APPT„_fc is presented. 
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3. On THE SPECTRUM OF LARGE WiSHART MATRICES AND RANDOM INDUCED STATES 

In the following, if (o„) and (6„) are some nonzero sequences, as usual, a„ ~ bn means that anjbn —>■ 1 
as n —>■ 00 , while o„ = o{bn) means that anjbn —>■ 0 as n —)■ oo; also, a„ ^ bn means that bnjan —)■ oo as 
n —>■ oo. 

Let us first recall the notion of Wishart ensemble of random matrices: 

Definition 3.1. A Wishart matrix of parameters (d, s) is a random d x d matrix W given hy W = GG*, 
where G € Aldxs(O) is a d x s matrix whose entries are i.i.d. complex Gaussian random variables of zero 
mean and unit variance. 


The following well-known result describes the behavior of the spectrum of a Wishart matrix of parameters 
(d, s) in the asymptotic regime d —>■ oo and s/d ^ c G (0, oo): 

Proposition 3.2. Let Ai > A 2 > • • ■ > Ad > 0 be the eigenvalues of a Wishart matrix of parameters (d, s). 
Then, in the asymptotic regime d —>■ 00 and s = 84 ^ cd for some constant c > 0, one has: 

(1) Almost surely, as d ^ 00 , the renormalized empirical eigenvalue distribution 

1 

Td = 

i^l 

converges weakly to the Marcenko-Pastur distribution 
TTc = max(l - c,0)(5o + \/4c- (x - 1 - c)^ 

(2) For any function jd = o(d), almost surely, as d ^ 00 , the rescaled eigenvalues Xi = d“^Ai have the 


following limits 


and 


Ad; Ad—1; . ■ . ; Ad—dQ-t-1 


0, ifc<l, 

(yc-l)2, ifol, 


Ai, A 2 ; ■ ■ •; Ajd -G be — {'/c + 1)^; 

(3) For every fixed fraction p G (0,1), almost surely, as d ^ 00 , the rescaled eigenvalue Aj^pdj converges 
to the (1 — p)-th quantile of the Marcenko-Pastur distribution: 


lim Aipdj — gi-p, 

d—y 00 


where qi-p is uniquely defined by 


diTc = 1 — p. 


Proof. The convergence in distribution stated in the first point is the classical result of Marcenko and Pastur 
[24]. The convergence of the extreme eigenvalues of Wishart matrices has been shown by Bai and Yin 
[6, Theorem 2], see also [5, Theorem 5.11]. The convergence of the eigenvalues in the bulk toward the 
corresponding quantile follows from the continuity of the distribution function of the limiting probability 
measure tTc (except at the eventual atom in 0), see also [27, Problem 2.4.19]. This also implies the convergence 
of the top (resp. bottom) o(d) eigenvalues towards the edges of the support of tTc- □ 


We consider next the asymptotic regimes where d <C s and, respectively, s <C d. 


Proposition 3.3. Consider a sequence of random matrices (Wg), where Wg G A4d, (C) is a Wishart matrix 
of parameters {dg,s), with dg = o(s). Then, 


Ve > 0, lim P 

s—>-oo 


El 

s 



< £ 


= 1 . 


( 1 ) 


Similarly, if (Wd) is a sequence of Wishart matrices of parameters (d, Sd), such that Sd = o{d), then 


lim P 


< e 

d—¥(X> 

d 



= 1 . 


( 2 ) 
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The first statement in the result above shows that in the regime d s, the eigenvalues of Wg/s converge 
to 1 as s —>■ oo, in probability (see [25] and the references therein). We recall the following result from [11], 
which shows that the fluctuations of these eigenvalues around 1 are semicircular. 


Proposition 3.4. Let Wd be a Wishart matrix of parameters {d^s) with s = Sd, and let 

Zd = 




be its centered and renormalized version. In the asymptotic regime 1 d ^ s (i.e. s/d ^ oo as d ^ oo), 
the random matrix Zd converges, in moments, to a standard semicircular distribution. Moreover, for any 
function jd = o{d), the top (resp. bottom) jd eigenvalues of Zd converge to 2 (resp. —2): almost surely, as 
d oo, 

XiiZd),...,XiiZd)^2 

Xi{Zd),...,Xi_^^^,{Zd)^-2. 

Proof. The convergence in moments has been shown in [11, Corollary 2.5], while the convergence of the 
operator norm of Zd to 2 has been shown in [11, Theorem 2.7]. The extension of the norm convergence to 
that of a o{d) number of eigenvalues is a classical argument in random matrix theory, see also the proof of 
Proposition 3.2 and [27, Problem 2.4.19]. □ 


We consider now the standard model of random induced quantum states. Let ^ be a random unit vector 
uniformly distributed on the unit sphere in 0 C'*. We denote by fj,d,s the distribution of the random 
quantum state 

p = Trc= {ipif*) e Dd, 

obtained after partial tracing over the ancilla space C'* the Haar-distributed random pure state tfip*, which 
is called the induced measure of parameters (d, s). Any random state p distributed according to the induced 
measure p,d,s is called a random induced state and we say that p comes from the induced ensemble with 
parameters {d,s). 

Random induced states are closely related to Wishart matrices. Indeed, it follows from [25, Lemma 1] 
that if W is a Wishart matrix of parameters (d, s), then p := is a random state with distribution pd,s- 
Therefore, results about Wishart matrices can be translated to random induced states, as we can see from 
the following result (for more details, we refer the reader to [2] or [25]): 

Proposition 3.5. IfW is a Wishart matrix of parameters {d,s), then for every e > 0, 


TtW 


ds 


- 1 


> e 


< C exp(—cdse^). 


for some C,c> 0. 


The main advantage of this approach is that the distribution of a Wishart matrix is much easier to deal 
with than the induced measure Pd,s- 


In the rest of the paper, we shall be interested in random bipartite induced quantum states: we shall 
assume that d = nk, such that we have a framework for studying the entanglement of the random states 
p ~ Pnk,s- Our techniques come from random matrix theory and are adapted to the study of large dimensional 
states {nk —)■ oo). We shall thus consider three asymptotical regimes: 

(1) The balanced regime: both n,k ^ oo; 

(2) The first unbalanced regime: n is fixed and k —>■ oo; 

(3) The second unbalanced regime: k is fixed and n —>■ oo. 
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4. Thresholds for RED in the balanced and the first unbalanced case 


We discuss in this section some questions left open in [18] regarding the thresholds for the set RED, 
in two asymptotic regimes. Indeed, in [18, Sections VII, VIII] it has been shown that for any size of the 
environment s which behaves like s ^ cnk for some positive constant c > 0, with overwhelming probability, 
quantum states p distributed along the induced measure with parameters (rik,s), will satisfy the reduction 
criterion, in the regimes where n grows and k is either fixed or it grows at the same speed as n (see [18, 
Theorem 7.2 and Theorem 8.2]). Hence, the thresholds for the set RED in these asymptotic regimes must 
be of smaller order than nk. We compute next the exact regimes and threshold values for these cases. 

Theorem 4.1. Consider a sequence of random density matrices (pn), where comes from the induced 
ensemble with parameters (nfc„,s„). In the balanced regime, where n —>■ oo, —>■ oo as n —>■ oo (not 

necessarily at the same speed) and Sn cn for some constant c > 0, we have (below, P = P„, as a function 
of n, denotes the probability distribution of pn): 

(1) If c < 1, then lim„_,.oo P[pn has positive reduction) = 0; 

(2) If c> 1, then lim„_>oo P[Pn has positive reduction) = 1. 

In other words, the threshold for the reduction criterion in the balanced regime is c = 1, on the scale Sn ^ cn. 


Proof. Instead of working with the induced measure for random quantum states, we shall use the simpler 
Wishart ensemble, since the reduction criterion is scale invariant: 

(id 0 R) > 0 ^ (id 0 R){W) > 0. 

To this end, consider a sequence (Wn) of Wishart matrices of parameters (nfc„,s„), and define 


Qn — RVl.n ^ ITn- 


Let us first assume that c > 1. Choose any e > 0 small enough such that c(l — e) > 1 + e and thus, for 
n large enough, s„(l — e) > n{l + e). From the partial trace property of Wishart matrices, it follows that 
the matrix WA,n follows a Wishart distribution of parameters (n, Since n = o{knSn), it follows by 

Proposition 3.3 (1) that 


lim P 

n—>-oo 


Wa.U C) hr. 


knSr, 


> (1 - e)In 


Similarly, using Proposition 3.3 (2), we have 


lim P 

n—>-oo 




= lim P 

71—^00 


Wn 

knSn 


< 


= 1. 


(3) 

n(I+e)^ 

^Tlkn 

= I. 

(4) 


We conclude in the case c > I by combining equations (3) and (4). 

Let us now move on to the case c < I, and show that, in this regime, with large probability as n —)■ oo, 
the random matrix Qn is not positive-semidefinite. We proceed in a similar fashion: using Proposition 3.3, 
we have 

lim P [WA,n C) 4„ < (1 + £)knSnInkJ = 1 (5) 

71—700 

and 

lim P[||W'„|| > (I - e)nfc„] = I. (6) 

71—700 

Choosing e > 0 small enough such that, for n large enough, (I — e)n > (I + e)s„, and using (5)-(6), we can 
conclude. □ 


Theorem 4.2. Consider the first unbalanced regime, where n and s are fixed integers, and k —>■ oo. Let 
(pk) be a sequence of quantum states, where pk comes from the induced ensemble with parameters (nk,s). 

(1) If s < n, then limfe_>ooP[pfe has positive reduction] = 0. 

(2) If s > n, then limfe_>oo P[pfe has positive reduction] = I. 

In other words, the threshold for the reduction criterion in the first unbalanced regime is s = n, on the scale 
of bounded s. Here, P = P^ is a function of k and denotes the probability distribution of p^. 
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Proof. The proof of this result is similar to that of Theorem 4.1. Working with Wishart matrices instead 
of random quantum states, the matrices Wk and WA,k have Wishart distributions of respective parameters 
{nk,s) and {n,ks). Assume first that s > n. From Proposition 3.3, it follows that 

lim P [WA,k 0 4 > (1 - e)kslnk] = 1 (7) 

k—^oo 

and 

lim P [IjVFfcll < (1 + £)nk] = 1. (8) 

fc—>-oo 

Choosing £ > 0 small enough such that (1 — s)s > (1 + s)n, and using (7)-(8), we can conclude. We leave 
the details of the case s < n to the reader. □ 

Note that the second unbalanced regime, where k is fixed and n —>■ oo, has been treated in [18, Proposition 
10.3]: the threshold, on the scale s„ ~ cnk, is given by 

{i + VkTTy 

k{k-i) ■ 

5. Thresholds for ARED 


Let us first recall some notations and results from [19]. For a vector x G K!j_, we denote by rkx the number 
of non-zero elements of x. 


Definition 5.1 {‘‘‘‘Hat operation” x i—> 5;). Given n,k > 2 and a vector x G with r < min(n, fc), we 
associate to x the pure quantum state i/; € C" 0 given by 

r 

■0 = ^ 0 /i, 

Z=1 


where and are fixed orthonormal families in C" and C^, respectively. We then define x to 

be the vector of eigenvalues of the reduction of the state ipif* G Dn,k, taken with multiplicities as 

in Corollary 3.3 in [19]: 


X := {xi,...,Xi,r]i,X2,...,X2,...,r]r-i,'- 


p, 0,..., 0 , Tjj,) s 


k—1 times 


k — 1 times 


k—1 times {n—r)k times 


where Xi > rji > Xi+i for i G [t — 1] and rjr = Vi A 0. The set {r]i}l^i \ equals the set of 

solutions rj gM. \ to the equation ~ ^ more details we refer the reader to Theorem 

3.1 in [19]). Moreover, if r = rka;, then —pr < (1 — 1/t) YZ.=i with equality if and only if Xi = Xj for all 
i,j G [r] (see Lemma 3.2 in [19]). 


A characterization of the set ARED„,fc has been given in Theorem 4.2 in [19] and states as follows: 


Proposition 5.2. We have 

ARED„_fe = {pG Dn,k\yx G A^in{n,k), > 0}, (9) 

where Aj; is the vector of eigenvalues of p ordered decreasingly and x^ is the increasingly ordered version of 
X that has been introduced in Definition 5.1. 


Let us make two general remarks about the set ARED„_fe, which provide upper and lower bounds for this 
set: 

{ A I At < (A: + l)At J C ARED„,fc C { A | rk A > (n - 2)fc + 2}. 

Lemma 5.3. For any probability vector A G ARED„ fc; we have rk A > (n — 2)k + 2. 

Proof. Although the statement follows from the inclusion ARED„_fc C LS 2 fc-i in [19, Theorem 8.1], we 
give here a direct proof. We assume n > 2 in order to avoid degenerate situations. Consider the vector 
a;=(l/2,l/2,0,0,...,0)e i,). Using the “hat operation” from Definition 5.1, we have 

f = ( 1/2,..^.,1/2 , ,-l/2)eM"'=. 

2k —1 times (n—2)k times 

If rkA < (n — 2)k + 1, then, obviously, {X^,x'^) = —At/2 < 0, and thus A ^ ARED„^fe. □ 
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Lemma 5.4. For any probability vector A G such that Aj < (fc -|- 1)A^^, we have A G ARED„_fe. 

Proof. For any x G ^min{n,k) of rank r = rkx, we have 


(A'^I 2;i (^Xnk + ■ ■ ■ + ^(n-l)fe-|-2) + ^l'^(n-l)fc+l 

+ X2 {^>^lri-l)k + -^ -^|n-2)fc+2) + V2>^\n-2)k+l 

+ Xr H -+ + ^'■^1 


>{k- l)(a;i H-h Xr)Xif, + (m + 

= ik-l)Xik + Vr{X\-Xik) 

>ik-l)Xi,-{l-l/r){X\-Xi,) 

>ik-l)Xi,-{l-l/k){Xi-Xi,) 


k-1 


ik + l ) Xi ,- X \ 


> 0 . 


+ Vr-l)X^f. + rjrXi 


By Proposition 5.2 it follows that A G ARED„ fe, which concludes the proof. 


□ 


5.1. Threshold for ARED in the balanced case. 


Theorem 5.5. Consider the balanced asymptotic regime, where n — )■ oo and fc„ — )■ oo as n ^ oo, and write 
s = Sn us a function of n. Let be a random induced state distributed according to the induced measure 
h-nk„,s,,- Almost surely, as n ^ oo and s„ ^ cnkn for some constant c > 0, 

(1) if c > 1, then Pn G ARED„ 

(2) if c < 1, then pn ^ ARED„_fc^. 


Proof. Let us start with the easier, second point: if c < 1 then, for n large enough, we have s„ < {n — 2 ) kn +‘ 2 , 
so, by Lemma 5.3, no eigenvalue vector A sampled from the induced measure with parameters {nkn, s„) will 
be an element of ARED„ (with non-zero probability). 

In the case c > 1, we show that the hypothesis of Lemma 5.4 is satisfied, almost surely as n —> oo. Indeed, 
using Proposition 3.2, we have the following almost sure limit: 


lim Sr, 

n^oo 




= Oc — 0 > 0, 


and thus, by Lemma 5.4, it follows that p„ G ARED„_/j^. 


□ 


The result above states that the threshold for ARED 

Sji ~ cnkn. 


in the balanced regime is c = 1, on the scale 


5.2. Threshold for ARED in the first unbalanced case. 

Theorem 5.6. Consider the first unbalanced asymptotic regime, where n is fixed and k ^ oo. Let pk be 
a random induced state distributed according to the induced measure Pnk,Sk ■ Almost surely, as k ^ oo and 
Sk ck for some constant c > 0, one has: 

(1) If c> n — 2, then pk G ARED„^fe; 

(2) If c <n — 2, then pk ^ ARED„^fe. 

In other words, the threshold for ARED in the first unbalanced regime is c = n — 2, on the scale Sk ck. 

Proof. Again, the second point follows from Lemma 5.3. Let us prove the first statement. To this end, we 
need to check that for all vectors x G A„, (A^^,i:^) >0. If rka: = I, the previous inequality is satisfied, since 
all the components of x are non-negative. Consider now a vector x with r := rka; > 2, and let e G (0,1) such 
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that c > n — 2 + e. Define tk '■= min(s/c, (n — l)fc). For k large enough, tk > {n — 2 + e)k. Then, we have 

+2^2(A|„_2)fc+2 + ^('n-2)fe+3 +-^ J 

>X\r]r+X2iek-l)Xj^ 

> x\r]r + X2{ek - l)Af^ 

>-At+n"^(efc - l)At^. (10) 


Note that the quantity appearing in the last step above is independent of a; € A„; we show next that this 
quantity is, almost surely as fc —>■ oo, converging to a positive limit. Define now, in the setting of Proposition 
3.2, for t > 0 


I ot, if t > 1 

at = \ 

\ai/t, ift<l 

to be the left-most positive element of the support of the free Poisson distribution ttj. Note that we assume 
t ^ I and thus d* > 0. 

Let us consider first the case c ^ n. Using the fact that, almost surely, as fc —> oo, SkX\ —> 6c/n 
SkX^^ —>■ ttc/n > 0, we conclude that the first negative term in the bound (10) vanishes, while the second 
term converges almost surely, as fc —>■ oo, to £j{nc)dc/n- 

Let us now treat the case c = n; this case requires special treatment because the left edge of the support 
of TTi is 0. Write, as before, for fc large enough. 


(Ap^,i^) > At77^ -b a;2(A|„_2)fc+2 + ^(n-2)fc+3 -^ 

> x\r]r + X2 L(fc — 1)/2J A|(„_ 3 / 2 )fe] 

>—X\ + n ^[(fc — l)/2jA|^^_3^2)fc]> 


where [•] denotes the ceiling function and [-J denotes the floor function. As before, we have that A^ —t 0, while 
■s/cA||-^_ 3 ^ 2 )fc] converges almost surely, as fc —>■ oo, to the 3/(2n)-th quantile of the free Poisson distribution 
TTi, which is positive (see Proposition 3.2). This concludes the proof. □ 


5.3. Threshold for ARED in the second unbalanced case. 


Theorem 5.7. Consider the second unbalanced asymptotic regime, where n —>■ oo, fc is fixed and Sn ~ cnk 
for some constant c > 0. Let {pn) be a sequenee of random states, where pn comes from the induced ensemble 
distributed according to the induced measure Pnk,s„ ■ Almost surely, when n —>■ oo and Sn ~ cnk, one has: 

(1) //c> (l + I -b l^/k + l)^, then pn G ARED„_fe; 

(2) //c < (l + I -b |Vfc + l)^, then ^ ARED„_fe. 


Proof. To prove the first statement we use again Lemma 5.4. Indeed, by Proposition 3.2, it follows that, 
almost surely, as n —)■ oo. 




(fc -b l)ac - be 


fc(v^)2 - 2(fc-b 2)v^-b fc > 0, 


if i/c > 1 -b I -b IVfc + 1 and thus G ARED„_fc. 

In the other case, for x = (1/fc, 1/fc,..., 1/fc) G A/j, we have 


d:= (l/fc,...,l/fc, ,l/fc-l)GK"^ 

k^ — l times {n—k)k times 


Thus, 

“ fc) + fc -^ a|„_j,)j,_^2) • 

Using again Proposition 3.2, it follows that, almost surely, as n —)■ oo, 

Sn (A^^, x^) [(fc -b l)ac - be ]. 

If c < 1, then Oc = 0 and hence the limit above is negative. On the other hand, if c > 1, then (fc -b l)ac — be = 
k{^/c)^ — 2(fc -b 2 )^/c + fc, which is also negative if c < (l -b | -b |-\/fc -b l) , and the proof is complete. □ 
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Remark 5.8. It is of interest to notice that when fc = 2, the threshold for ARED in the second unbalanced 
case is c = 7 + 4\/3 and it coincides with the one obtained in the same regime for APPT [11]. This result is 
natural since the two criteria are equivalent when the second subsystem is a qubit (see [19] and the references 
therein). 


6. Thresholds for LSp 

In [19], the authors introduce for each p G [nk], the set of eigenvalue vectors for which the largest eigenvalue 
is less or equal than the sum of the p smallest: 

LSp := {A G A„fc : X\ < + X^,,_p_^_^ H-h (11) 

It is worth to mention that the set LSp is of particular interest because it sets bounds for (the more compli¬ 
cated set) ARED„_fc- Indeed, accordingly to [19], it holds that for n, fc > 2, 

LSfc C ARED„,fc C LS2fc-i. (12) 

Note that the set LSp depends only on the product d = nk, and not on the particular values of n and fc. 
We compute now the threshold for the sequence of sets in three different cases: pd = p is a, hxed 

function, 1 Pd = o{d) and pd = \ td \, for some Hxed fraction t G (0,1). 

Theorem 6.1. For d = nk, let pd be a random induced state distributed according to the induced measure 
Pd,s- Then, in the asymptotic regime d —>■ oo and s = Sd cd for c > 0, one has: 

• For every fixed integer p > 2, almost surely, as d ^ oo, 

(1) if c> {1 + then pd G LSp/ 

( 2 ) ifc<(lF^^^Y,thenpd^TSp. 

• For every function pd such that 1 <C Pd = o{d), almost surely, as d ^ oo, 

(1) if c> 1, then p G LSp^; 

(2) if c < 1, then p ^ LSp^. 

• For fixed t G (0,1) and pd = [td\, almost surely, as d ^ oo, 

(1) ifol- t, then p G LSp^j ; 

(2) */ c < 1 - t, then p ^ LSp^j. 

Proof. We start with the case when p is fixed. For every fixed p > 2, the inequality from (11) becomes 
asymptotically be < pOc, where the constants Oc and be are defined in Proposition 3.2. In the case when 
c < 1, Uc = 0, so the inequality above cannot be satisfied. In the other case, when c > 1, the inequality is 
equivalent to c > (1 -|- Y, and the conclusion of the theorem follows. 

We move now to the second case where pd —t oo as d —)■ oo, with pd = o(d). If c > 1 then, almost surely, 
by Proposition 3.2, 


SdAf —>■ be 

SdX'^, ..., SdA^_p^_,_^ —>■ Oc > 0. 

Using Pd oo, we obtain Sd{X\ -I- • • ■ -b A^_p^_,_^) —>■ oo, finishing the proof of the first point. On the other 
hand, if c < 1, the limiting measure tt^ has an atom at zero of mass 1 — c. Since the number of eigenvalues 
appearing on the right hand side of the inequality defining LSp_j is pd = o(d), these eigenvalues will be zero 
(for d large enough), while SdX\ be > 0, proving the claim. 

Let us now consider the final case, where the parameter pd behaves like td, as d — >■ oo (t G (0,1) is 
fixed). Assume first that c > 1 — t. Notice that, since A^ converges to zero, it is enough to show that 
Ai + --- + Al [tdj-i-i converges to some positive limit. To show this fact, let us consider two sub-cases. First, 
if c > 1 > 1 — t, then, as d —)■ oo. 


A^ H-1- > L^dJA^ it/c)ae > 0. 
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Otherwise, we have 1 > c > 1 — t. Take £ > 0, small enough such that c — 3£ > 1 — t. Then, we have, for d 
large enough, 

-^ ^d-itd\+i = -^ >'i-itd\+i 

- '^t(c-e)dj -^ >^d-[td]+l 

> [edj . 

The right hand side of the above expression converges almost surely, as d —)■ oo, to a positive constant {e/c 
times the (1 — c + £)-th quantile of tTc, which is positive), proving the first point. 

We consider now the second point, when c < 1 — t < 1. Since, almost surely, rkp^ < Sd ~ cd, we have 
= • • • = = 0, while SdX\ ^ 6c > 0, showing that pd i LSp^j. □ 


7. Thresholds for GER 


In [19], the following set of probability vectors (we put r = min(n, k)) was introduced: 

{ r —1 r —1 'j 

AS : ^Ai'<2At, + ^Ai,_,,|, (13) 

in connection to Hildebrand’s characterization of APPT states. Indeed, it was show in [19, Theorem 7.2] that 
GER„ fc C APPT„ fc; the proof consists in applying Gershgorin’s circle theorem to show that Hildebrand’s 
conditions from [13] are satisfied. 

We compute next the thresholds for the set GER„ fc. Note that the definition of the set is symmetric in 
n and k, so we shall assume, without loss of generality, that n > k. 


Theorem 7.1. Consider the balanced asymptotic regime, where n —>■ oo and —>■ oo as n —>■ oo. Let pn 
he a random induced state distributed according to the induced measure Pnk„,s„- -Almost surely, as n —> oo 
and Sn ~ cnk^ for some constant c > 0, one has: 

(i) If c> A, then S GER„,fe„ ; 

(ii) If c < A, then ^ GER„,fe„. 


Proof. We consider the eigenvalues of the corresponding Wishart matrix with parameters {nkn,cnk^). Note 
that the ratio of parameters is ckf —)■ oo as n —>■ oo, so we can apply Proposition 3.4. We shall consider the 
rescaled eigenvalues 




From Proposition 3.3 we know that, for all i, Ai —>■ 1 in probability, while Proposition 3.4 implies that, for 
any function = o{nkn), almost surely as n ^ oo. 


V Z — 1,2,..., jk; '\/c^n(Aj 1) 2 

i — l;2,...,jVi; '\/c^n(Anfc,i + l —i 1 ) ^ 2. 


Let us first assume c > 4. The inequality (13) for the Xi is implied by the following equivalent inequalities 
(recall that min(n, kn) = kn)'- 

{h-n l)Al A (^n T I)Xnk„ 

kn{Xl Xnkjf) A Ai + Xnkn- 

The left-hand side of the inequality above converges, almost surely, as n —)■ oo, to Aj^/c, while the right-hand 
side converges to 2; since c > 4, the conclusion follows. 

In the case c < 4, we write 

Ai -l- • • • -l- Afc„_i > {kn — l)Afc^_i and 
2A„fc„ -|- Xnkn-l + • • • + Xnk„-k„ + l C: {kn + l)A„fe„_fe„ + l. 

Using Proposition 3.4 with j„ = kn = o{nkn), we can conclude as before. □ 
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Theorem 7.2. Consider the unbalanced asymptotic regime, where k is fixed and n —>■ oo. Let pn be a 
random induced state distributed according to the induced measure Pnk,s„ ■ Almost surely, as n ^ oo and 
Sn cnk for c > 0, one has: 

(i) If c> {k + — ifi, then pn € GER„^fc; 

(ii) If c < {k + y/k'^ — 1)^, then pn ^ GER„^fc. 

Proof. As usual, since the definition of the set GER is scale-invariant, we shall consider a Wishart matrix of 
parameters (rik,cnk). The inequality from (13) reads now asymptotically (recall that k is fixed) 

{k - l)hc <{k + l)ac, 

where the constants a^, are defined in Proposition 3.2. If c < 1, Oc = 0, so the inequality above cannot be 
satisfied; we assume thus c > 1. The inequality is then easily seen to be equivalent to c > (fc -|- \/fc^ — 1)^, 
and the proof is complete. □ 

Remark 7.3. The thresholds for GER are the same as the thresholds for APPT given in [11, Theorem 4.1] 
and [11, Theorem 4.2]. This result strengthens the claim that GER is a very good approximation to APPT. 
Similar results to support this claim are [19, Remark 7.4] and [19, Proposition 8.2]. 

8. Threshold for SEPBALL 

One of the earlier results about the geometry of the set of separable states is a very surprising one: the 
largest Euclidean ball, centered at Id/d, which is contained in the set of d-dimensional quantum states, 
contains only separable states [12] (to make a sense of separability, we consider an arbitrary decomposition 
= C" 0 C^). To be more precise, the set 

SEPBALL„,fc = |p e 1 Tr(p2) < 

is a subset of the set of separable states. Moreover, since the set SEPBALL is defined in terms of the 
trace of the square of the density matrix, it is invariant under global unitary conjugations, so we have that 
SEPBALL„,fe C ASEP„,fe. 

It turns out that the set SEPBALL is much smaller than the other sets studied in this work. In [19, 
Proposition 8.2] it has been shown that the largest eigenvalue of elements in SEPBALL„^fe is smaller than the 
corresponding quantity for other sets, such as APPT„_fc, ARED„_fe, or GER„^fe. The behavior of thresholds 
is also different for SEPBALL„ ^ than for the other sets: the size of the “environment” Sd scales like the 
sguare of the total size of the system d = nk. In the result below, since SEPBALL„_fe depends only on the 
product d = nk, we simply write SEPBALL^; = SEPBALL„ fe. 

Theorem 8.1. Consider the asymptotic regime, where the total dimension d = nk of the system grows (the 
way n and k grow is not relevant). Let pd be a random induced state distributed according to the induced 
measure Pd,sd ■ Almost surely, as d ^ oo and Sd cd^ for c > 0, one has: 

(i) If c> 1, then pd € SEPBALL^; 

(ii) If c < 1, then pd ^ SEPBALL^. 

Proof. Replacing pd by Wd/PrWd, where Wd is a Wishart matrix of parameters (d, Sd), the inequality in the 
definition of the set SEPBALL^ reads 

(d-l)Tr(VEj) < [TrWdf. 

In our setting, d = o{sd), so we can apply Proposition 3.4. Replacing Wd by Sdid + y/dsdZd, the above 
inequality simplifies to 

dss{d - l)Tr(Zj) < ds^ + 2sd\/d^Tr0d -|- dsd(TrZd)^. 

Since Zd converges almost surely to a semicircular distribution, we have that d~^TT:Zd —t 0 and d“^Tr(Zj) —>■ 
1. Hence, the previous inequality becomes, after replacing Sd by cd^ and keeping only the dominating terms 
in d —)■ oo. 


finishing the proof. 


cd® < c^d^. 


□ 
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Balanced regime 

Unbalanced regime 


n, fc —>■ 00 

m = min(n, fc) hxed, max(n, fc) —>■ 00 

SEP 

< s < log^ n [4, n=k] 

mnk ^ s ^ mnfclog^(nfc) [4] 

PPT 

s ~ cnk 

s ~ cnk 

c = 4 [2, n = fc] 

c — 2 -|- 2 -\^ 

1-7U [8] 

RLN 

s ~ cnk 

s hxed 

c = (8/37r)^ [3, n = fc] 

s = w?' [3] 

RED 

s ^ cn 

For TO = n, s is hxed 

For m = k, s = cnk 

c = 1 

s = n 

[18] 


Table 1. Thresholds for separability vs. entanglement and entanglement criteria 


9. Conclusions and open questions 

In this final section we gather results about the thresholds for different entanglement criteria considered 
in the literature and also for some related sets in the balanced and unbalanced asymptotic regimes. 

In Table 1 we review the thresholds for separability (SEP) vs. entanglement and also for some well- 
known entanglement criteria: positive partial transpose (PPT) criterion, realignment (RLN) criterion and 
reduction (RED) criterion. In [3] the authors showed that, in the balanced case (n = fc —>■ oo), the threshold 
for the realignment criterion is c = (S/Stt)^ « 0.72, on the scale s ~ cnfc, which means that the realignment 
criterion is asymptotically weaker than the PPT criterion (from a volume perspective; from a set-inclusion 
perspective, the two criteria are not comparable). In the unbalanced regime it is shown that the threshold 
is s = 771^, on the scale of bounded s; here, m = min(n, fc). On the other hand, the corresponding one for 
the PPT criterion is unbounded with respect to max(n, fc), s ^ cmax(n,/ c)to. Eor the reduction criterion, 
the unbalanced case splits into two different cases depending on the parameter which tends to infinity (the 
dimension of the subsystem to which the reduction map is applied is important, for more comments see [19]). 
When the dimension k of the second subsystem is larger, then the threshold is s = n (n is the dimension of 
the first subsystem), on the scale of bounded s. In the other case (n —oo and k is hxed), the threshold for 
the reduction criterion is c = on the scale s ~ cnk, which is smaller than the corresponding one 

for the PPT criterion, c = 2 -|- 1 — which follows from the paper of Banica and Nechita [8]. Moreover, 

for k = 2 the two values are the same. This is natural since the reduction criterion is in general weaker than 
the PPT criterion, the two criteria being equivalent for k = 2. In [23], upper bounds for the threshold of 
/c-extendible states have been obtained, but we do not discuss these results here, since the fc-extendibility 
criterion does not enter the framework of this work. 

In Table 2 we gather the thresholds for the set of absolutely PPT states (APPT) from [11] and for GER 
and ARED. Eor k = 2, the threshold value computed for ARED reads c = 7-|-4\/3; this value coincides with 
the one obtained for APPT, which is in agreement with the fact that APPT„ 2 = ARED„ 2 - The thresholds 
for GER are the same as the thresholds for APPT, which shows that GER is a very good approximation to 
APPT (see also [19]). Since the sets SEPBALL and LSp depend only on the product d = nk, it is sufficient 
to consider only one asymptotic regime (d —>■ 00 ), and thus we consider a separate table (Table 3). 

Let us finish this work with a list of open questions: 

(1) Find a description of the set ARLN of quantum states satisfying the absolute version of the realign¬ 
ment criterion; following [1], this is a superset of APPT. Gompute the thresholds for the set ARLN 
in different asymptotic regimes. 
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Balanced regime 

Unbalanced regime 


n. A: —>■ oo 

m = min(n, k) fixed, max(n, fc) —>■ cx) 

APPT 

s ~ cmin(n, k)'^nk 

s 

~ cnk 

c = 4 [11] 

c = (to + Vto^ — 1) [11] 

GER 

s ~ cmin(n, k)'^nk 

s 

~ cnk 

c = 4 

c = (to 

+ ^/m? — 1)~ 

ARED 

s ~ cnk 

For m = n, s ^ ck 

For m = k, s = cnk 

c = 1 

c = n — 2 

c = (1 + I + + 1) 


Table 2. Thresholds for related sets 



Asymptotic regime d = nk ^ oo 

SEPBALL 

s ~ cd'^ 

c = 1 

LSp 

s ^ cd 

p>2 fixed 

1 <C p = o{d) 

p= [td\, te (0,1) 

VP-i) 

c = 1 

c = 1 — t 


Table 3. Thresholds for SEPBALL and LSp 


(2) Give a simple description (or tight bounds) for the set ASEP. Using this description, compute the 
thresholds for the set, in different asymptotic regimes. The values of these thresholds, compared to 
those for APPT, could invalidate the conjecture [1] that ASEP = APPT. 
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